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Reassessing nuclear matter incompressibility and its density dependence 
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Experimental giant monopole resonance energies are now known to constrain nuclear incompress¬ 
ibility of symmetric nuclear matter K and its density slope M at a particular value of sub-saturation 
density, the crossing density pc. Consistent with these constraints, we propose a reasonable way to 
construct a plausible equation of state of symmetric nuclear matter in a broad density region around 
the saturation density po- Help of two additional empirical inputs, the value of po and that of the 
energy per nucleon e(po) are needed. The value of K{po) comes out to be 211.9 ± 24.5 MeV. 
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I. INTRODUCTION 

The nuclear incompressibility parameter Kq defined 
for symmetric nuclear matter (SNM) at saturation den¬ 
sity Po stands out as an irreducible element of physical 
reality. It has an umbilical association with the isoscalar 
giant monopole resonances (ISGMR) in microscopic nu¬ 
clei; it also underlies in a proper understanding of super¬ 
nova explosion in the cosmic domain [l|. From careful 
microscopic analysis of ISGMR energies with suitably 
constructed energy density functional (EDF) £{p) in a 
non relativistic framework as applicable to finite and in¬ 
finite nuclear systems, its value had initially been fixed 
at Kq — 210 ± SOMeV [ 3 , 0 . In microscopic relativistic 
approaches on the other hand, a higher value of Kq ^ 260 
MeV was obtained . After several revisions from differ¬ 
ent corners, however, its value settled to Kq ~ 230 ± 20 
MeV [i^. It gives good agreement with the experimen¬ 
tally determined centroids of ISGMR, in particular, for 
^°®Pb , ®°Zr and ^^"^Sm nuclei, calculated both with non- 
relativistic [1, 0 and relativistic 0 energy density func¬ 
tionals. The near-settled problem was, however, left open 
with the apparent incompatibility of the said value of 
Kq with the recent ISGMR data for Sn and Gd-isotopes 
0, [l0l - [T7j| . These nuclei showed remarkable softness to¬ 
wards compression, the ISGMR data appeared explained 
best with Kq ^200 MeV 

A plausible explanation was recently put forward by 
Khan et.al [l8l| for the apparent discrepancy. It is ar¬ 
gued that there may not be an unique relation between 
the value of Kq associated with an effective force and 
the monopole energy of a nucleus predicted by the force 
M- The region between the center and the surface of 
the nucleus is the most sensitive towards displaying the 
compression as manifested in the ISGMR. The ISGMR 
centroid Eg is related to the integral of incompressibility 
(f K{p)dp) over the whole density range [2^. As a result, 
a larger value of K{po) for a given EDF can be compen¬ 
sated by lower values of K (p) at sub-saturation densities 
so as to predict a similar value of ISGMR energy in nuclei. 
It is seen that the incompressibility K {p) calculated with 
a multitude of energy density functionals when plotted 


against density cross close to a single density point [l8l| . 
this universality possibly arising from the constraints en¬ 
coded in the EDF from empirical nuclear observables. 
This crossing density pc[= (0.71 ± 0.005)po] HH seems 
more relevant as an indicator for the ISGMR centroid. 
Because of the incompressibility integral, the centroid 
seems more intimately correlated to the derivative of 
the compression modulus (defined as M = 3pK'{p) ) 
at the crossing density rather than to Kq. The value of 
Kc{= K{pc)) is seen to be ~ 35 ±4 MeV [2l|. From var¬ 
ious functionals, the calculated values of Mc{= M{pc)) 
are found to be linearly correlated with the correspond¬ 
ingly calculated values of ISGMR centroids for ^°®Pb and 
also for ^^°Sn. From the known experimental ISGMR 
data for these nuclei, a value of Me ~1050 ±100 MeV 
is then obtained, revised from an earlier estimate of 
1100 ± 70 MeV [l^. Using a further assumption of a 
linear correlation between Kq and Eg calculated from 
different EDF, a value for Kq ~230 MeV with an un¬ 
certainty of ~40 MeV is reported, the uncertainty being 
inferred from the spread of Kq values obtained with the 
different functionals used. 

The universality of the crossing point pc and the val¬ 
ues of Kc and Me can be readily acknowledged; Me is 
seen to be well correlated to Eg- The Pearson correla¬ 
tion coefficient r of Me with Eg for ^^'’Sn is 0.80 
and is 0.94 for ^°®Pb. However, assumption of a lin¬ 
ear correlation between Kq and Eg may not be justi¬ 
fied, they seem to be very weakly correlated (r=0.67 for 
^^°Sn and 0.79 for ^°®Pb) [2lj. The inferred value of 
incompressibility around saturation may then be called 
into question. One can see that a linear Taylor ex¬ 
pansion Ko{pq) = K{pc) ± (po - Pc)K'{pe) yields for 
Ko — 185 ± 14.3 MeV, noting that K'{pe) = Me/{3pc). 

The absence of a strong linear correlation between Kq 
and Eg calculated from different effective forces prompts 
one to think that Ke and Me alone are not sufficient to 
yield the correct value of Kq. Further empirical informa¬ 
tion is possibly needed to arrive at that. In this paper, 
we show that with given values of only Ke and Me along 
with some time-tested values of empirical nuclear con¬ 
stants, it is possible to address to a proper assessment of 
the value of incompressibility K and its density depen- 
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dence. The empirical constants are the saturation density 
Po, taken as 0.155 ± 0.008 fm“^ for SNM and the energy 
per nucleon at that density e(po}, taken as —16.0±0.1 
MeV [ 2 ^, [ 2 J. An acceptable value of the effective nu¬ 
cleon mass m*jm, which lies in the range m*/m ~0.8 
±0.2 [ 2 ^ at saturation density is also used. 

This paper is structured as follows. In Sec. II, we in¬ 
troduce the theoretical elements to calculate the nuclear 
equation of state from and pc with the aid of empirical 
inputs mentioned. Results and discussions are presented 
in Sec. III. Sec. IV contains the concluding remarks. 


II. THEORETICAL EDIFICE 

We keep the discussions pertinent for SNM at any den¬ 
sity p at zero temperature (T = 0). The chemical poten¬ 
tial of a nucleon is given by the single-particle energy at 
the Fermi surface, 


, = = ^ + U (1) 

where ppip) is the Fermi momentum and U{p) the single¬ 
particle potential. Assuming the nucleonic interaction to 
be momentum and density dependent, the single-particle 
potential separates into three parts 

U = Vo+p%Vi + V2. (2) 


The last term V 2 is the rearrangement potential that 
arises only for density-dependent interactions, and the 
second is the momentum-dependent term that defines the 
effective mass m*, 


— = — + pIVi 

2m* 2m ^ 


so that 


A = - + 2^i- 

m m 

The energy per nucleon at density p is given by, 
p^ 1 2 1 




m , p 

—) < 
m 




(3) 


(4) 


(5) 


The density dependence of the effective mass [l^ can 
be cast as ^ = 1 ± /cp, the rearrangement potential can 
be written in the form V 2 = ap°‘. This is the form that 
emerges for finite range density-dependent forces [ 2 ^ in 
a non relativistic framework or for Skyrme interactions. 
The quantities a, a and k are numbers. If ^(po) is 
chosen, k is known. 
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At p = Po, P = 0, then from Eq. (7), writing for ^ = 


&p2/3 with6=(M^, 


eo = e(po) = ^Po^^[I-2A:po]-apE^. (8) 

Since P = p^|^, from Eq. (7) again we get. 

At Po, this yields (since Ao = 

\aap^ ± ^hkp^^ “ ^ 

Furthermore, Eq. (9) gives 

K{p) = 9^ = 26p2/3 - I65fcp3/3 
op 

d^P 

-9a(Q; ± l)ap“ ± 9p-—(11) 

ap‘‘ 

Defining M = 3p^ = 27p^^, this leads, at p = pc to 

9a(a ± l)ap“ ± l&hkpl/^ - {2hp]l^ + ^ - K^) = 0(12) 

Since fc is a given entity and pc and (Mc/3 — Kc) are 
known, eqs. (8) and (12) can be solved for a and a, 
eq. (10) then gives the value of the nuclear incompress¬ 
ibility Kq. Once Kq is obtained, Mo{= M(po)) is 
evaluated from eq. (12) by choosing po for pc- Then 
Qo = 27pg|^|pp is also known from Mq = 12Aro ± Qo- 
The structure of eq. (9) shows that the pressure and its 
first derivative are interrelated. One can then get higher 
density derivatives of P or of energy e recursively from 
eq. (9) as is evident from eq. (11). For the present, we 
show that 


From Gibbs-Duhem relation, 

P 

p = e± —, 
P 


( 6 ) 


where P is the pressure. Keeping this in mind, from 
Eqs. (1),(5) and (6), we get 


, \ Pp 2to, P 

e(p) = tt;—[3-^ - E 2 H-, 

IUto m* p 


(7) 


9p^-^ = 9a^(a ± l)ap“ ^ ± ^bkp^/'^ - Uyp (13) 
Op'^ 3 3 


Since 


we find 


where we have put < p^ >= |pp. 


Q‘ip d'^e „ d^e 


d^P 1 

OPo'^^lpo = ^ATo ± 2Qo ± g-^o- 


(15) 
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where we have defined Nq = Slpg^lpp. From eq. (13) 
and (15), knowing Kq and Qo, Nq can be calculated. 
Similarly, one can calculate the fifth density derivative of 
energy {Rq = 243pg|^|pJ by exploiting eqs. (13) and 
(14) from 

d^P 8 1 

9Po'^^lpo = 4(5o + giVo + —i?o- (16) 

These help to find the density variation of the energy and 
also of the incompressibility, as is seen, 

e(p) = e(po) + 

2 D 

where e = () (counting terms only up to e® is seen to 
be a very good approximation in the density range of ^ 
po/4 < p < 2.0po, we retain terms up to them). Eqs. (7) 
and (17) give 

^ = e(po) + + iQoe^ + ^Noe^ 

p z D z4 

+ Y^-Roe®-^P^^^[l-2A:p] + ap“- (18) 
and eq. (9) gives 

K{p) = 9^ = 18[— - + \bkp^/^ + iaap“](19) 

dp p 15 3 2 

We have thus the equation of state (EOS) of symmetric 
nuclear matter in a reasonably spread-out density domain 
around the saturation density. 

The incompressibility K at any density p can be cal¬ 
culated directly from eq. (19) or it may be calculated in 
terms of K{pc) and its higher density derivatives as 

K{p) = K(p,) + (p-p,)7^'(p,) + i^_^iF"(p,) 

D 

The different derivatives can be calculated from eq. (19). 
With given values of po,eo, ^(po), and pc, one notes 
that the solutions for a and a do not depend separately 
on Kc and Me, but on {Mc/3 — Kc). 

III. RESULTS AND DISCUSSIONS 

The values of the empirical constants po, eg and ^ 
needed for our calculation have already been mentioned. 
As for the crossing density, we choose pc = 0.110 ±0.0008 
fm“^. With given inputs of and Kc, it should be 
noted that the output values for Me and Ke may come 
out to be different, but {Mc/S — Kc) remains invariant. 



FIG. 1: (Color online) The sensitivity of the incompressibility 
Ko at the saturation density (po) on the values of the incom¬ 
pressibility Ke (green dash-dotted line), its density slope Me 
(red dashed line), the crossing density pe (blue dotted line) 
and the value of po (black full line). The abscissa extends 
from -1 to +1. These end points refer to the scaled lower and 
upper limits of Kc, Me, Pc and po, respectively (see text). 


With inputs Mc=1050 MeV and Kc =35 MeV, the out¬ 
put Me and Kc are found to be 1051.8 MeV and 35.46 
MeV, respectively. Since they are very close to the input 
values, they were not tinkered with for exact matching 
of the output and input values. The value of incompress¬ 
ibility at Po turns out to be iCg=211.9±24.5 MeV either 
from eq. (19) or eq. (20). We note that in eq. (20), at sat¬ 
uration, the value of the second term on the right hand 
side is 143.3 MeV, the third term is 35.9 MeV, the fourth 
term is —3.2 MeV, the fifth term (not shown in eq. (20)) 
is 0.55 MeV and so on,which adds up to ^ 211.9 MeV. 

The uncertainty in an observable X (like K, M etc) is 
calculated from where Ayi are the 

uncertainties in the empirically known entities yi. The 
sensitivity of Kq on these entities that influence the in¬ 
compressibility most is displayed in FiglTJ The abscissa 
is scaled such that 0 refers to the central value of these 
entities Me, Kc, Pc and pg; ±1 refer to the extrema of 
their domain (±100 MeV, ±5 MeV, ±0.005pg and ±0.008 
fm“^ from the central values of the entities, respectively). 
The value of Kq is seen to be very sensitive with changes 
in either Me or pg when all other input entities are kept 
fixed. Its sensitivity to Kc or pc is weak; on ^ or to the 
energy per nucleon eg, it is rather insensitive. The near¬ 
insensitivity of incompressibility to the effective mass is 
observed for Skyrme density functionals also. From the 
data base for these functionals as tabulated by Dutra et. 
al [ 1 ^, the correlation coefficient between Kq and m* is 
calculated to be only ^ —0.2. 
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FIG. 2: (Color online) 

The incompressibility and its different density derivative as 
defined in the text plotted as a function of Me- 

The near-perfect linear correlatioir of Kq with Me as 
seen in FiglT] is very startling. From Eq. (20), one 
may expect that the second and higher order deriva¬ 
tives of K{pe) would destroy this correlation. However, 
we find that both K" and K'" are also linearly corre¬ 
lated with Me and thus K{po) retains its linear corre¬ 
lation with Me- This is displayed in FiglU where we 
define Ki = {po - pe)K'{pe), K 2 = K"{pe) and 

_ (po-^pe) x"'{pe). The weak correlation between Kq 
and Me that can be inferred from the calculated correla¬ 
tion structure of {Me — Ea) and {Kq — Eq) in refs. MM 
possibly results from the use of different EDFs in getting 
the various relevant observables. 

Figures [3] and 0] display the functional dependence of 
the nuclear EOS on density. The panels (a) and (b) in 
FigElshow the energy per nucleon and the pressure, the 
ones in Fig0] show the incompressibility and its density 
derivative M, respectively. As one sees, the uncertainty 
in energy and pressure grows as one moves away from 
the saturation density, similarly the uncertainty in in¬ 
compressibility or its density derivative increases with 
distance from the crossing density. 

IV. CONCLUSIONS 

To sum up, we have made a modest attempt to re¬ 
assess the value of K{po) consistent with the new-found 
constraint on the incompressibility K{pe) and its den¬ 
sity slope M{pe) at a particular value of density at sub¬ 
saturation, the crossing density pc- We have relied on 
some empirically well-known values of nuclear constants. 
We have further made the assumption of linear density 



' I I I I — I _ ,_ I_,_□ 

0 0.5 1 1.5 2 


P/Po 


FIG. 3: (Golor online) The nuclear EOS as a function of den¬ 
sity. The panels (a) and (b) show the energy per nucleon and 
pressure, respectively in a selected range around the satura¬ 
tion density. 



FIG. 4: (Golor online) The nuclear EOS as a function of den¬ 
sity. The panels (a) and (b) show the incompressibility and its 
density derivative M, respectively in a selected range around 
the saturation density. 


dependence of the effective mass and the power law de¬ 
pendence of the rearrangement potential which happens 
to be generally true for non-relativistic momentum and 
density dependent interactions. In relativistic models, 
the density dependence of the effective mass may not be 
linear (^ . The rearrangement potential appears explic¬ 
itly there only in the case of density dependent meson 
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exchange models [2^. 

The value of incompressibility K{pq) turns out to be 
211.9 ± 24.5 MeV. This is somewhat lower than the cur¬ 
rent value in vogue, Kq ~ 230 ± 20 MeV. From recursive 
relations, our method allows also estimates of higher den¬ 
sity derivatives of energy or of pressure and thus helps 
in constructing the nuclear EoS e{p) at and around the 


saturation density. 

The authors gratefully acknowledge the assistance of 
Tanuja Agrawal in the preparation of the manuscript. 
One of the authors (JND) acknowledges support from 
the Department of Science & Technology, Government of 
India. 


[1] G. E. Brown, Phys. Rep. 163, 167 (1988). 

[2] J. P. Blaizot, Phys. Rep. 64, 171 (1980). 

[3] M. Farine, J. M. Pearson, and F. Tondeur, Nucl. Phys. 
A615, 135 (1997). 

[41 D. Vretenar, T. Niksic, and P. Ring, Phys. Rev. C 68, 
024310 (2003). 

[5] B. G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett 
95, 122501 (2005). 

[61 P. Avogadro and C. A. Bertulani, Phys. Rev. G 88, 
044319 (2013). 

[71 T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C 78, 
034318 (2008). 

[81 B. K. Agrawal, S. Shlomo, and V. K. An, Phys. Rev. C 
72, 014310 (2005). 

[9] J. Kvasil, D. Bozhik, A. Repko, P.-G. Reinhard, 
V. Nesterenko, and W. Kleinig, arXiv:1412.1997 [nucl- 
th] (2014). 

[10] Y. W. Lui, D. H. Youngblood, Y. Tokimoto, H. L. Clark, 
and B. John, Phys. Rev. C 70, 014307 (2004). 

[11] D. H. Youngblood, Y.-W. Lui, B. John, Y. Tokimoto, 
H. L. Clark, and X. Chen, Phys. Rev. C 69, 054312 
(2004). 

[12] T. Li and et. al,, Phys. Rev. Lett. 99, 162503 (2007). 

[13] J. Piekarewicz, Phys. Rev. C 76, 031301(R) (2007). 

[14] U. Garg and et al. T. Li, Nucl. Phys. A788, 36c (2007). 

[15] V. Tselyaev, J. Speth, S. Krewald, E. Litvinova, 
S. Kamerdzhiev, N. Lyutorovich, A. Avdeenkov, and 
F. Grmmer, Phys. Rev. C 79, 034309 (2009). 

[16] T. Li and et. al„ Phys. Rev. C 81, 034309 (2010). 


[17] D. Patel and et. al,, Phys. Lett. B 718, 447 (2012). 

[18] E. Khan, J. Margueron, and 1. Vidana, Phys. Rev. Lett. 
109, 092501 (2012). 

[19] G. Colo, N. V. Giai, J. Meyer, K. Bennaceur, and 
P. Bonche, Phys. Rev. C 70, 024307 (2004). 

[20] E. Khan, J. Margueron, G. Colo, K. Hagino, and 
H. Sagawa, Phys. Rev. C 82, 024322 (2010). 

[21] E. Khan and J. Margueron, Phys. Rev. C 88, 034319 
(2013). 

[22] S. Brandt, Statistical and Computational Methods in 
Data Analysis (Springer, New York, 3rd English edition, 
1997). 

[23] B. K. Agrawal, J. N. De, and S. K. Samaddar, Phys. Rev. 
Lett. 109, 262501 (2012). 

[24] N. Alam, B. K. Agrawal, J. N. De, and S. K. Samaddar, 
Phys. Rev. C 90, 054317 (2014). 

[25] M. Dutra, O. Lourenco, J. S. SaMartins, A. Delfino, J. R. 
Stone, and P. D. Stevenson, Phys. Rev. C 85, 035201 
( 2012 ). 

[26] D. Bandyopadhyay, C. Samanta, S. K. Samaddar, and 
J. N. De, Nuclear Physics A 511, 1 (1990). 

[27] A. Bohr and B. R. Mottelson, Nuclear Structure, vol. Vol. 
I (Benjamin, New York, 1969). 

[28] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 
(1986). 

[29] S. Typel and H. H. Wolter, Nucl. Phys. A656, 331 
(1999). 



